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,  20.  Abstract  (continued) 

are  described  and  compared  to  data  from  the  recovered  experimental  model.  The 
tunnel  closure  predictions  and  the  predicted  locations  of  block  motion  are 
in  reasonable  agreement  with  some  experimental  data. 

Other  aspects  of  the  study  included: 

•  development  and  application  of  an  analytical  technique  for 
predicting  joint  activation  conditions  (i.e.,  block  motions) 
for  specified  joint  orientations  and  locations  near  the 
tunnel  surface, 

•  calculations  and  comparisons  of  tunnel  response  to  dynamic 
and  quasi-static  stress  loads,  and 

•  investigation  of  the  effects  of  joint  friction  on  tunnel 
damage. 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 


1 . 1  BACKGROUND 

All  rock  masses  contain  large  numbers  of  planar  discontinu¬ 
ities  or  joints.  These  joints  often  react  to  an  applied  stress 
quite  differently  than  the  surrounding  material  and  there  is  both 
experimental  (Ref.  1)  and  theoretical  (Refs.  2,  3)  evidence  to 
suggest  that  sliding  motion  across  these  joints  should  be  one  of 
the  important  considerations  involved  in  designing  underground 
structures  which  are  to  survive  nuclear  and/or  conventional 
explosions. 

In  the  recent  DIABLO  HAWK  underground  nuclear  test,  SRI 
International  fielded  several  laboratory-scale  structures  in 
jointed  rock  simulant  (Ref.  4).  Figure  1  shows  the  general  model 
designs.  The  jointed  rock  models  are  18"  and  30"  in  diameter, 
with  3"  and  5"  diameter  tunnels,  respectively.  The  SRI  small- 
scale  jointed  structures  experiment  provides  an  opportunity  for 
examining  the  influence  of  highly  idealized  jointed  arrays  on 
tunnel  closure  and  damage  characteristics. 

There  are  both  continuum  and  explicit  approaches  for  pre¬ 
dicting  tunnel  closure  and  damage  in  jointed  media.  In  the 
continuum  approach,  the  overall  mechanical  properties  of  the 
field  arc  synthesized  from  the  properties  of  the  intact  rock  and 
the  properties  and  arranqcment  of  the  joints;  for  example,  stiff¬ 
nesses  and  strengths  used  in  the  continuum  model  are  lower  than 
the  intact  rock  yiropcrtics.  In  the  explicit  approach,  the  prop¬ 
erties  and  arrangement  of  individual  blocks  and  joints  are 
individually  and  explicitly  modeled  and  used  to  calculate  tunnel 
response  to  loading.  Doth  of  these  treatments  have  important 
limitations:  The  rationale  for  choosing  a  set  of  degraded 


a 


I 


propur  Lies  in  the  continuum  approach  is  usually  hard  to  justify, 
and  this  approach  precludes  discontinuities  in  displacements 
across  joints,  which  may  be  a  critical  omission.  The  properties 
of  joints  in  the  explicit  approach  are  hard  to  measure,  and  the 
arrangement  of  joints  in  a  real  field  around  a  real  tunnel  can 
only  bo  estimated. 


l.i  OBJbu'rr\'H  AND  APPROACH 


I'n..,  'bnctive  of  this  study  is  to  examine  the  influence  of 
ti'.e  )oinLS  on  tunnel  closure  and  damage  using  numerical  techniques 
whicii  explicitly  treat  the  idealized  joints  in  the  SRI  structure. 

The  app'roach.  involves  three  related  aspects: 

i.  Joiiit  .'vcti  vation  . 

An  analytical  closed  form  solution  is  derived  for  tiie 
initiation  or  activation  of  relative  motion  across  a 
joint.  Tnc  Joint  activation  depends  on  (a)  the  rela¬ 
tive  orient^ition  of  the  joint  and  tlic  tunnel  surface 
intersected  by  tile  joint,  (b)  tlie  magnitude  and  propa- 
aation  direction  of  the  stress  wave  encompassing  tiic 
tunnel,  u-ki  (c)  tiio  siiear  stress  on  the  joint  whicii  mu.st 
excecu  a  critical  value  wiiich  depends  on  the  normal 
stress  across  the  joint  as  well  as  on  surface  roucjlincss. 

[n  Section  d,  tile  joint  activation  solution  is  Jerivi,wi 
aiii,i  ap'pd  it'd  to  the  SRI  small-scale  models  to  estimatt' 

(a)  tiie  iireas  on  tiie  tunnel  surface  wdiere  joint 
actiViition  is  likely  to  occur,  (b)  the  stress  levels 
wiiii’ti  are  necessary  to  activate  joints  of  various 
or  1  (.nita  1 1  ons  ,  and  (c)  the  sensitivity  of  joint 
rj'Ct  1  v,:i  L  ion  to  material  [iro;  >er  t  ios  . 
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2 .  Dynamic  Tunnel  Kesponse  Cases . 


In  the  SRI  small-scale  models  with  tunnel  diameters  of 
3" ,  the  transit  time  across  this  characteristic  dimen¬ 
sion  is  only  22  nsec  in  the  16A  rock  simulant  material. 
Tims,  since  the  rise  times  to  peak  stresses  of  0.8  kbar 
in  the  DIABLO  HAWK  event  are  5-10  msec,  these  small- 
scale  structures  are  actually  loaded  (and  unloaded)  in 
a  quasi-static  fashion.  Full-scale  tunnels  may  be 
loaded  more  dynamically,  therefore  the  effects  of 
dynamic  loadinq  on  small-scale  tunnels  buried  in  jointed 
media  are  analyzed  using  six  plane  strain  computer  code 
simulations  with  dynamic  loading  times  of  100-400  i:sec. 
These  cases  also  vary  the  assumed  joint  and  tensile 
strenqtii  proijorties  of  the  rock  simulant.  Tlie  WAVE-L 
Lagranqiian  code  is  used;  the  basic  formulation  is 
similar  to  the  HEMP  code  (Ref.  a)  .  Tiie  numerical 
solutions  and  material  properties  are  di..scnbod  in 
Section  J. 

J.  Ojioji  i  qd  t'T  t_i  c'  'luiinel  !<..'Sponse  Prediction  For  a  Sma  1 1-Sca  le 
S^r'Jotu_re  Loaded  to  C) .  8  k_b_a_r  Pur ina  the_  DIABLO  liAWK 
!•', veil  t . 

A  numericai  imulat  luii  i^rcuiction  of  tunnel  closure  ami 
a.i:i:a  :e  is  made  and  compareii  with  a  fleldeFl  small-scale 
L I  I'eee!  (  I  de  ;nis  a  3"  tunnel  dianK'ter  with 

.naai-.d  iti  -  -aMuticn  tecnnique  and  thti  c  a  i  culaL  i  ona  i - 
■  i  i;','_nt  il  ccm',;.  ai  iseiis  are  dt_'scribed  in  S'l-ction  4. 


1.3  SUMMARY  AND  CONCLUSIONS 
i.3.1  Joint  Activation 

Althoucjl:  the  response  of  a  tunnel  in  a  jointed  medium  must 

be  solved  numerically,  the  conditions  that  will  initiate  joint 

slipping  can  be  analyzed  in  closed  form  if  wo  assume  an  elastic 

media  with  joints,  and  a  far  field  stress  state  of  uniaxial 

strain  with  a  stress  loading  of  o, 

load 

Figure  2  shows  the  loading  stress  required  to  activate 

a  joint  oriented  at  ,t  =  0,  ‘45°,  and  90°  as  a  function  of  the 
angle  (^^)  wnere  the  joint  intersects  the  tunnel  surface.  Thus 
a  joint  oriented  at  i  =  45°  will  require  a  loading  stress  of 
.06  kbar  (  900  psi)  to  become  activated  if  this  joint  intersects 
tile  tunnel  surface  at  3  =  100°  (see  the  following  sketch  and 
Figure  2 )  . 


Ik-liavior  for  !'■  +  180° 

N 
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Note  that  for  each  joint  orientation  (')>  the  sliaded  reqi(3n 

on  Figure  2  shows  the  angles  (.)  on  the  tunnel  surface  where  joint 

activation  is  i)Ossib2e  as  a  function  of  loading  stress  ( 

Note  that  synuiiotry  implies  equivalent  tunnel  behavior  at  angles 

of  .  +  180".  Joints  will  not  slip  outside  of  these  shaded  regions. 

The  tunnel  angles,  ,  associated  with  absolute  minimums  in  load- 

inn  stress  (>  ,  )  necessarv  for  joint  activation  are  indicated 

’  load 

in  Table  1.  However,  note  that  the  activation  stress  curves 
(F'ic;uro  2)  are  relatively  fl.at  bottomed;  and  therefore,  joint 
activation  can  occur  over  a  fairly  wide  range  in  tunnel  angles 
with  only  a  small  increase  in  loading- 

Taiiie  1.  Tunnel  AiiLiles  (,  )  for  Minimum  Loading 

Stress  Lor  Various  Joint  Orientations  (  <) . 


1  .Joint 

'  Orientation 

'funnel 

IJ 

Angies  lor  Minimum 
oaiiing  Stress* 

'  1 

( 

^  ciL-’ijroc 

0  i 

1  ,  ■' 

ut-'gree 

,  ,  load 

Kbar 

psi 

07,  113 

.068 

9  8G 

i  -43 

103 

.0  59 

8  3b 

;  0  ; 

i 

41,  139 

.ltd  2 

330 

, .. 

a  >  1 

_  i 

7o,  103 

.039 

S3b 

*l':c-;)r>.n;inLaL  i  ve  constitutive  properties  are  usi'd 
(das'.'  I'l)  :  Joint  Coliesion  (io)=-U'l  kbar.  Friction 
Angli.'  (:)  il".  Material  Poisson's  Hatio  ( \“ )  =  •  2 13  . 

Figuia;  >  stiows  Liie  areas  on  L l\e  tunnel  surface  whe're  joint 
ict,  ivition  (liu’  a  single  joint)  coulil  occur  for  ,  =  .3  kbar. 
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All  the  joint  orientations  in  the  SRI  experiment  are  shown  for  a 
representative  set  of  joint  and  eartli  media  properties.  (Varia¬ 
tions  in  the  joint  and  earth  media  constitutive  properties  are 
illustrated  in  Section  2.) 

1.3.2  Dynamic  Response  Cases 

Six  WAVE-L  numerical  simulations  (Table  3)  were  performed 
involvinq  100-400  iisec  stress  rise  times  to  a  peak  loading  stress 
of  O.S  kbar  on  the  SRI  3"  diameter  models.  The  wave  transit  time 
across  the  3"  tunnel  opening  is  22  iisec,  therefore,  the  dynamic 
tut\nol  response  characteristics  could  be  investigated. 

I'Ue  li,'tu;iLe  stress  history  near  the  tunnel  crown  is  strongly 
a<-'i  etuleii  t  '.'ll  the  stress  .vave  rise  time  for  times  corresponding 

to  a  ti-'W  ir.iiiSLt  times  across  the  tunnel  diameter.  A  stress 

w.ii'e  M.-a,'  t  i::ie  el  100  nsec  ('-5  transit  times)  caused  a  peak  ten- 
lie  ..a  I  e:...;  ;  .ilr.iest.  1  j  0  Durs.  In  the  guasi-static  loading  case 
wi'.nea  ii.i  iiaina,  no  tensile  stress  developed  near  the  crown. 

. '  .1.  ;  ■■  •;  i.n.in  :  et  tile  material  near  the  tunnel  surface  is 

■  ilee  1.  i  I'l  ;  !  a>.  I.  1  a  I  1  ve  d  i  spl  acemen  ts  *  across  the  joints. 

Vd*'  :;  .  ir;  ■  ■.•ai  a...-re,ss  )oiuts,  tlui  plates  of  rock  separated 

;  ;  !  .e  .  i:;'.;  ,m:i  .e.M  as  unsuppoi' teti  btiams  due  to  the  j.'resence 
e!  ise  ..uiiii.e.;  i.aUK.l.  "n.  resu  1  I  i  iia  iji.vim  can  K\ui  to 

‘  e  ■ : . :  I  1  .a  .  .  1 : .  1  M  .  ic  t  u  re  . 

1  .  1  .  i  i.i  a  a  s  1  -  t  a  t  1  e  Tunnel  Response  Rrediclion 

liu.-  ;  a.  1 . .  1  -  s  t  a  t  1  e  solution  tr-chnigue  is  ilescribOLl  in  Section  4. 
a  '  r  1 , -a  1  sirailat  ion  ptedii.'t  ilUks  usinq  this  teclinigue  were  made 
!  a  .1  !  1  e  1  lod  small-scale  (}"  d  i  .atK' t r  tunnel)  SRI  model.  The 

*  1  no  1  e  1  a  t  1  \'e  1 1  s;  1 ,  i  oiuiu  ui  t.  or  slij'ping  .icross  a  joint,  is  measured 

h;  tiiv-  cuitont  distance  betweiui  t.wo  points  which  were  initially 
toU'Ctiina  oacn  otiioj  on  oippos  i  te  laces  of  the  joint. 


calculation  predictions  were  made  prior  to  excavation  and  examin¬ 
ation  of  the  model. 

Figure  4  shows  a  comparison  of  the  predicted  tunnel  closure 
(as  a  t  of  the  diameter)  versus  tunnel  angle,  .  The  agreement 
is  reasonable  except  near  the  crown- invert  (i  =  0°  and  180°) . 

The  maximum  predicted  closure  is  1.2°  at  c-  =  30°.  The  experi¬ 
mental  data  show  a  maximum  1.6%  closure  at  the  crown-invert 
compared  to  a  0.5°6  predicted  value  at  the  crown-invert.  This 
difference  may  be  due  to  the  moderate  sliding  on  the  joints  and 
associated  tensile  failure  and  relaxation  which  occurs  in  the 
calculation,  but  not  in  the  experiment. 

In  the  post-shot  experimental  model,  there  was  deformation 
of  a  thin  (.003  inch)  stainless  steel  tunnel  liner  which  surrounds 
the  foam  filler.  This  liner  may  have  inhibited  sliding  along  the 
joints.  The  observed  liner  deformation  is  associated  with  the 
joints  which  are  predicted  to  slide  the  most;  thus,  some  "block 
motion"  appears  to  have  occurred. 

It  is  concluded  that  tunnel  closure  in  a  well-characterized 
jointed  media  could  bo  predicted  using  existing  computer  codC' 
techniques  involving  explicit  descriptions  of  joint  locations, 
orientations,  and  friction  properties.  Also,  because  of  the 
important  role  joint  (and  fault)  motion  plays  in  determining 
damage  to  tunnel  liners  and  neaiby  oartli  material,  it  is  felt 
that  continuum  approaches  to  predicting  tunnel  closure  and  damage 
might  not  bo  as  useful. 
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SECTION  2 


JOINT  ACTIVATION  ANALYSIS 


The  necessary  conditions  for  a  joint  or  block  to  begin  slip¬ 
ping  depend  on  the  relative  , :  L  of  the  joint,  the  tunnel 
surface,  and  the  stress  wave,  as  indicated  in  the  following 
sketch : 


\ 

\ 


I  is  the  joint  Oiientn- 
tion  angle 

is  Che  tunnel  angle 
(angle  where  joint 
inLersi.'1-ts  tunnel) 

is  Che  angle  between 
die  joint  and  die 
normal  to  the  tunnel 
sur  i.iee  ,  .  =  1  ! 


Joint  sii[jjMc)e  will  occur  when  the  shear  stress  (i)  on  the  joint 
cxceei^ls  a  critical  value,  i  .  This  value  t-ici’cnds  on  the  normal 
stress  (  jj)  across  the  joint  as  well  as  on  the  surface  roughness, 
i’ri.i^r  to  loint  activation,  the  earth  media  around  the  tunnel  is 
assumeei  to  respond  elastically  and  statically  to  a  uniaxial  strain 
louii  from  a  (' luasi-static )  pL^inar  wave. 


1  (. 


Thu  key  ciiaracteristics  of  interest,  with  respect  to  Joint 
are  : 


•  the  areas  on  tlie  tunnel  surface  where  joint 
activation  is  likely  to  occur, 

•  the  stress  levels  which  are  necessary  to 
activate  joints  of  various  orientations, 

•  the  sensitivity  of  joint  activation  to  material 
properties . 


Tne  determination  of  these  characteristics  involves  the 
foilowin>.j  elements: 


1.  The  •  tnnoontini  a  ' ac  field  on  the  uniined 

tunnel  surface  (this  stress  field  is  shown  in 
Figure  5  as  a  ratio  witii  tlie  loading  stress, 

i.e.,  . 

2.  The  i (i.e.,  we  assume 

I  =  i  +  (tan  .‘.)c  where  t  and  g  are 
max  o  n  o  ’ 

experimentally  determined  parameters) . 


J  . 


The 


o n  , 


:  'on  as  a  function  of  joint  orientation 
(  0 ,  position  on  the  tunnel  (2) ,  and  material 
projjer ties  . 


The  activation  stress  load  relationship  (derived  in 
Appendix  A)  IS 


wliere 


a 

load 


( 2  ;  ^cos  ;. )  /  ( sin  ( 2  V- ;  )  -  sin.*.) 


2  -  11  +  2cos  (2;-.)  ] 

1 1- V) 


(1) 


(sec  sketcn  on  previous  patje) 

i  and  ;  are  joint  constitutive  parameters 
o 

'  is  the  Poisson's  Ratio  for  the  earth  material 
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Note  that  Equation  1  shows  that  u^,  »  and  v  are  the  con¬ 

stitutive  properties  of  interest  for  joint  activation.  Table  2 
lists  tht'  numerical  values  of  those  properties  which  wore  used 
to  calculate  joint  activation  relationships. 


Table  2 .  Joint  and  Earth  Media  Constitutive 
Properties 


Joint 


Case  Cohesion*,  [ 

c 

kbar  psi 


Earth  Media 


Friction  Angle*  Poisson's  Ratio 


1  (degrees) 


*Tne  joint  cohesion  ( and  friction  angle  (:)  are  related  to 

the  unconfined  compressive  strength  of  the  earth  media  by 

2  :  cos  ; 
o _ 

a  I 

I  -  sin: 

Figure  2  and  Fiemres  ii  to  9  show  the  variations  in 

^  load 

versus  tunnel  angle  ( ;• )  for  the  joint  and  earth  media  property 
variations  indicated  in  Table  2.  Figure  3  and  Figures  10  Oe  li 
show  the  corresponding  regions  on  the  tunnel  surface  where  joint 
activation  is  possible  for  e,  ,  =  .5  kbar. 


loint.  Activation  Stress  Loads  tor  Various  .kiiut  i.)r  iontat  ions  (;)  as  a  iiiiution  ol 
iuiiiifl  Anglo  (■  )  lor  Constitutive  i’ropcrty  Caso  1’2  {:  =.\)S,  ;  -  ii"  ,  .  =  .2j;. 


(OEG)  BETA  tOEG) 


BETA  (DEG) 


Comparisons  of  the  results  from  the  constitutive  property 
cases  PI  to  P5  leads  to  the  following  conclusions: 

1.  Increasing  (or  decreasing)  the  cohesion, 
increases  (or  decreases)  the  minimun  loading 
stress  for  joint  activation  while  leaving 
unchanged  che  tunnel  angles  (B)  where  joint 
activation  is  most  liliely. 

2.  Reducing  the  fiction  angle  to  =  0,  thereby 

simulating  a  von  Mises  limiting  shear  strength 
of  causes  an  expansion  of  the  possible  joint 

activation  region. 

3.  Varying  Poisson's  Ratio  for  the  earth  media  from 
V  =  .25  to  V  =  .40  causes  an  expansion  of  the 
possible  joint  activation  region. 
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SECTION  3 


DYNAMIC  RESPONSE  CASES 

In  the  DIABLO  HAWK  experiment,  the  stress  wave  time  history 
involves  rise  times  to  a  peak  stress  of  0.8  kbars  in  b-10  msec. 
In  Section  4,  quasi-static  techniques  are  used  to  treat  this  rise 


Time  Ailer  Wave  Arrival  .U  O.S  khar  Stress  Rain;e 


time,  which  is  (|uite  icnij  with  respect  to  tunned  resi'onse.  In 
Section  3.1,  material  properties  for  all  tiie  cases  are  descnbcel. 
In  Seettion  3.2,  six  numerical  ca  1  e:u  1  a  t  ion  s  involvimi  much  sliorter 
rise  times  of  between  100  and  400  esec  ar<.'  LiescribcLi  and  comiJared 


3.1  MATERIAL  PROPERTIES 


IbA  Rock  Simulant 

In  tile  small-scale  experiments,  the  jointed  media  is  16A 
rock  simulant.  The  elastic-plastic  proiicrties  are  modeled  with 
the  following  parameters: 

Bulk  Modulus  K  =  .131  Mbar 

Shear  .Modulus  G  =  .086  Mbar 

Poisson's  Ratio  v  =  .23 

Density  =  2.0  gm/cm-’ 

Yield  Strength  Y  =  .22  kbar  +  1.2P,  where  P 

IS  pressure  .  A  Prandtl-Reuss 
non-associated  flow  rule  is 
assumed . 

Tensile  failure  is  modeleci  in  most  of  the  solutions  using  a 
tensile  failure  model  which  permit.s  oriented  cracks  to  develop 
(Ref.  b).  A  300  psi  =  34.5  bar  tensile  strength  was  measured 
in  a  simple  tensile  test  (Ref.  7)  on  a  16A  rock  simulant  sample. 
Most  of  our  solutions  use  tliis  value;  in  one  of  the  cases,  the 
tensile  strength  is  reduced  to  6.9  bars  (100  psi)  as  a  sensitivity 
study . 

la  two  of  the  cases,  frictional  forces  along  the  sliding 
joints  are  included.  In  these  situations  the  joint  constitutive 
model  is  specified  by: 

=  +  =  :  +  ( tan ;  )  ' 

ma.\  o  n  o  n 

w  here 

-  ma.ximum  allowable  sliear  stress  on  joint 

ma  X 
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=  normal  stress  on  ioint 
n 

r  =  the  ioint  cohesion 
o 

U  =  (tan.f)  is  the  coefficient  of  friction  and  ij) 
is  the  angle  of  friction 

The  experimental  values  (Ref.  7)  are 

T  =10  psi 
o 

-  30° 

The  Drucker-Prager  plastic  yield  surface  assumed  for  the  16A 
rock  simulant  material  does  not  include  shear  strength  softening 
(i.c.,  reduction  of  the  yield  strength)  due  to  plastic  shear 
deformation  or  tensile  cracking.  This  may  be  significant  in  this 
application  because  the  uniaxial  strain  loadintj  path  for  this 
material  is  close  to  the  yield  surface.  Hierefore ,  small  changes 
in  the  failure  surface  may  lead  to  large  changes  in  the  plastic 
response  ot  the  lOA  rock  simulant. 


Foam  Filler 

The  tunnel  is  filled  with 
with  the  following  properties: 

Bulk  Modulus 

Sliear  Modulus 

Density 

Von  Mises  Yield 
Streir  :th 


a  low  density  elastic-plastic  foam 

K  =  1.75  kbar 
G  =  2.5  kbar 
=  .05  gm/cm’ 

Y  =  500  psi  =  0.34  kbar 
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3.2  CASE  COMPARISOMS  FOR  THE  0.8  KBAR  DYNAMICALLY  LOADED  TUNNELS 


Table  3  lists  tlie  joint  geometry,  tensile  limit,  rise  time, 
and  material  model  for  the  six  dynamic  loading  cases.  In  all 
cases  the  tunnel  is  filled  with  low  density  foam  with  properties 
described  in  Section  3.1. 

The  loading  is  a  plane  wave  of  0.8  kbar  magnitude.  The 
stress  increases  linearly  from  0  to  0.8  kbar  over  a  characteris¬ 
tic  rise  time  i^,  {100  psec  in  most  of  these  cases)  .  A  constant 

pressure  of  0.8  bar  is  maintained  after  t  .  The  far-field  con- 
dition  of  the  pulse  (i.e.,  prior  to  interaction  with  the  tunnel) 
is  uniaxial  strain. 

3.2.1  Effect  of  Loading  Rise  Time 

In  the  first  three  cases,  the  rock  media  is  modeled  as  an 
elastic-plastic  media  witb  no  joints  or  tensile  failure.  Cases 
1,  2,  and  3  differ  in  their  loading  rise  times  (i^  =  100,  200, 
and  400  psec) .  Cases  4,  5,  and  6  have  joints  which  are  perpe: 
dicular  to  the  wave  vector  as  indicated  in  Table  3;  in  these 
cases,  the  joint  properties  and  tensile  strength  of  the  media 
are  varied. 

Figure  14  siiows  particle  velocity  (upper  part  of  figure) 
and  principal  stress  (lower  part)  fields  at  62  iisec  for  Case  1. 
Note  that  the  velocity  and  stress  fields  iiave  a  plane  of  symmetry 
wliicli  sei’arates  tile  two  fields.  The  scale  bars  for  the  velocity 
vectors  aiiLi  principal  stresses  are  stiown  on  tiiis  f  i<.;ure .  The 
following  si\etcii  liefines  the  plottin<i  convention  for  the  princiiial 
stress : 


! '  I  i  :i '  i  : '  I  1  I  1'' 


I' I  i  111  - !  I'.i  1  Si  I  1 

1 1 ' 1 1  ■  i  I'll 

1  ' 'wn-.' .  1 1'li  1  I  eipiiiu'ii  t  s 

,  I  r  >  r  I  >ri  I'l'i'  ill') 

1  1  ni|'  Ti-  S'.  1  I'll 


i  I  i  I  S  S  r  i  II 1  i  p.i  I 

I  t-s  . 
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Table  3.  The  Dynamic  Loading  Cases 


'I 


a 


Vf loci  tv  FiolJ 


">0  tt/soi' 


I'ortirlc  Vi’locitv  anil  Principal  Stress  l  icliis  at  bP  .,scc  lor  tlic 
i.  1  as  t  i  c-i' 1  as  t  i  c  .'-ioilia  Model  '.('nso  I,  O.S  khar  I.oadiiu;  in  11)0  ;  sc 


no  loading  wave  approachon  the  Lunno'i  from  L  ho  loft.  IV/  G2  viscc, 
no  Lanaiol  is  rospoiKiinp  very  non- sy iiu'io t r icai  1  y  .  Noto  that  ono 
rinoinai  ntronn  is  noarly  zoro  noat  thi.’  tui\noi  boundary,  ans-l  tlio 
amsontial  stross  is  iainjosL  near  tho  s.  ■!  i  n,;  in  no  ot  tiio'  tunnol. 

i  ioati-.'  14  otiows  tiuit  one  of  tho  [irinoipal  stinnsscs  is 
ii  ;n.tiy  L^-nisiio  near  Luo  tunnol  crown  (loft  si^io  of  tunnol)  . 
i.ruio'  Lj  snov.s  flio  tonsiio  strouss  hisLijiy  at,  t.iio  crown  for 
.iso.-  i,  1,  a'ui  J.  'I'ho  tousilo  ro‘s;H>ns<.,’  nocir  Lno  crc>wn  is  clearly 
■..':’.,.wn  t  >)n  tno  stress  wave  riSv.’  t  i  n,'.  for  tit.os  oi  a  fou  hundroii 
icr-.)LWConas  ana  ivass;  a  sounii  wav>.'  i-rossa.-s  th.o  cinj r ac  to r  i  s Li c 
'.nnua  a  1  1  i_:r  i.:'  >ibout  2.1  nsoc.  In.  to-.  a.i  s  .1  -  s  t  i  t  i  c  solution 

itnouL  annt..  activation  (si..'o  .so^'f  i,i!;  4)  ,  t  noi'..  was  no  t^unsilo 
t.  rose  ai  too  on 'wn. 

. .  t  t  ' c  t  .a  w’ ri  .  ‘ n. . o  It.,’  .a.  t  t  ■  n  t  n  ^u n.  . K' •  i  ! 1 1  .f i lo a  i'  ti  t I'o n a  t ii 

,’nt  ■  .'..hi  i .  - . ,  c. 'U'O  i:.'.  .ca;..a'ity  1 1 1  ..iiL.,.w  tiW:  formation  o; 

■.a. ill.  cr.tcns,  to  vioic...  tio  tin.'  wia;n  o !'  t  iio  clocks,  an.,!  to 

..oast  stro^oos  .i..).ic>'n’  i  :  On.  ci.u'ks  usi-i.;  clastic  s t n.' s s- s t ra i n 
o  L  .1 1  i '.;o  sn  !  os  .  .'•stt-'i  a  c>.  i  1  nts  ^'I'lcuoa,  ■  tonsil..'  stross  across 

:o-  era..':.  IS  oaimittoo,  .;ia  no  .nioar  .Or^-ssos  arc  p'..aoa,  L  t>..'d  on 
co.ica  o  i.iu'  cr  icK  IS  opn-n  (t.iio  v.  i,;t;i  .d  tin.'  cracn  is  con- 
!.nao,;si'.'  ol' ;  i  1 1  o  roa  )  .  dun-n-  stress  .ia  i  ind.io'iit  s  ,una  im.'.irtant 
'  ‘  ^n  i i  So.'  ,i  ..'lai.'.a  !  I't-.'ay.io  n  L  I  y  .iltt.‘rs  t  Ik  lo.o.il  str.i.'-.s  !  in  S'^icii 

w  1  y  a.,  t  '  oanini.’o  it.-;  own  .•in.wtn;  t.nus,  ur.'  tv-al  i:-;tic  niothoti 
:  :  in  ao.-t;:i  i  ■■i.n,',-.  ari'-..;tii  must  considi.'t  liio  alt’,  iw-d  stress  f  l  o  Ki  . 

i.'.i  -'.'li  Cato-nil  .-a  I  i  IS  ill.o'.-.o.l  tv!  tl.voU'i  up  to  two  cracks. 

!  I  0.1  rn  ■.'i  a'K  t.'i’.is,  tiia’  cv.'l  1  IS  '..‘v  >  n . 1  vit.!  I' o  Vi  tv'  i.'o  comp  I  O  t  a' 1 
a  1 o  t'-  .ui'.i  no  t'-'nsi.;>n  .it  .i  1  i  is  Liu' lain  f  to  r  ai  lowi.'vi.  Pno 
■  o..a,i;  t  tn,.'.  i  tmv’n  I  '  ;  t.v,-n;;  i  b.-  t.iiJairo  ari.i  tiio  :  lottiuia  aion  Va'n  t  i  on  s 

.....  t 1  in.iu;  t'.'  t  h-o  v't  aa'Ka'd 


anill.s  in  tno  v'oni.' 'Uf  a  L  i  ona  1  ari  i 


criticcii  tensile  failure  criteri(3n 
in  the  indicated  principal  stress 
direction ,  >  . 

Platie  whicii  exin-riencc-s  tensile 
fai lure . 


Crack  e,|jt_.ns  and  relieves  the 
tensile  stress  such  that  nc'  lorces 
or  Lrvictions.  exist  cai  the  crack 
surfaces;  i.e.,  Se'ro  normal  and 
shear  stress  on  the  crack  sur'f^ice. 


In  the  directions  perpendicular 
to  the  crack  surfaci.',  the  material 
c>in  maintain  ti-’iisile  stresses 
until  tlie  critical  tensile  failuri.: 
criterion  is  exceeded  aiul  a  maw 
crack  is  formed. 


P  1  ot.  1.  nu  Con  veil  t.  ion  : 

le  crack  in  x-y  plane.  Tlie  line  is  plotteti  in 
di  recti. Ill  cd  tlie  civick  (normal  to  tlui  principal 
ss  ixiusina  tile  cr.K'k)  . 


el-plaiie  i-iae’s;  1  ^  ^  _  c.iused  th.e  I'l'ack. 


'-plane  and  x-y  plane  ct  ick. 


x-y  plane  ctacks;  only  (lu'-et-pl ima  tension  is 
i  t  t  ed  . 


'etely  tia.-tuie.i  -  ne  teiisii.n  permit  te.,!. 


he 


I  I  a  at  i-;,n  t  ..1  lensili-  iia.k-  in  WAVI-1.  (.'.In. 


In  Cases  4  and  b  the 


/  of  the  solution  to  varia¬ 


tions  in  the  tensile  streiiijth  of  the  rock  medium  is  examined. 

In  Case  4,  the  IbA  rock  simulant  has  a  tensile  strenytti  of 
b.9  bars  (=  iOQ  [jsi)  ;  in  Case  5,  a  stremjth  of  34,5  bars 
(=  500  psi)  is  used.  'Die  effect  on  the  tensile  crack  formation 
IS  sliovvn  in  I'iqures  17  and  18  which  show  these  two  cases  at  two 
different  times. 

Case  fa  impiements  all  the  material  and  joint  properties  used 
in  thd  t|uasi-static  DIABLO  HAWK  simulation  (Section  4),  includinq 
friction  aloiu;  the  joints.  Tlic  presence  of  the  friction  is  the 
only  uifference  between  Case  5  and  fa.  'I'lie  final  crack  pattern  in 
Case  6  is  compared  to  Case  5  in  Fi<jure  19.  /•  /' '  '  •,  '• 

Tensile  crackiiui  of  tiie  material  near  the  tunnel  surface  is 
closely  rel^ited  to  tiie  relative  slide  dis[j  lacements  across  the 
joints.  Wiien  sii;.)piiuj  occurs  across  pOnts,  Lhi.-  pdates  o'  roc;; 
se[jarated  o'.'  tne  joints  can  act  ^is  unsupported  beams  vju<-  to  tl'io 
presence  of  tiie  unitned  tunnel.  'I'he  resultiiu;  beam  bt,‘n..iinq  leai,ls 
to  tensile  stresses  <,i!Ki  fracture. 
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SECTION  4 


DIABLO  HAWK  SMALL-SCALE  STRUCTURE  PREDICTION  CASE 

This  numerical  simulation  was  intended  to  match  a  fielded 
small-scale  SRI  model.  The  calculation  results  were  obtained 
■  to  the  excavation  of  the  model.  Tlie  particular  model  is 

designated  as  S-Ul-2-5  by  SRI. 

4.1  PROBLEM  CONEIOUIUM’ION 

The  tunnel  diameter  is  .  Joint  spacing  is  nominally 
(0.47")  .  The  boundaries  of  the  numerical  tirid  are  at  9.45"  from 
tne  tunnel  center  (u.i  times  the  tuiniel  radius) ,  which  is  suf¬ 
ficient  to  minimise  edue  effects  in  the  recjion  of  interest.  I’he 
joint  tieonietry  is  sliown  witii  tiie  initial  com',iuta t i ona  1  erii.i  in 
E'luure  20. 

Material  properties  of  tlie  IbA  rock  simulant  and  foam  filler 
.ire  di.'scribed  in  Section  3.1.  The  frictional  jiroin'rties  of  the 
joints  are  also  discussed  in  that  section. 

The  tunnel  and  jointed  media  are  subjected  to  a  planar  stress 
wave  of  0.8  kbnr  peak  am(>lituae.  The  rise  timi'  is  5-iil  r.sec 
lied  ih'ca’.  InU  m.Sec.  Tnest.-  tim.js  nri.'  very  lonn  comp.na'a  t  e, 
tile  22  ..sec  .sound  wave  transit  t.inic  across  tile  tunnel  diameter 
in  the  lOA  ruiterial.  Tile  actual  tunnel,  therefore,  is  in  a 
" ' ;uas  1- s  ta  t  ic "  st^ite. 

4.2  OLABI-STATIC  SOLUTION  TECHNICLJ: 

'[’lie  juasi-static  solution  is  fcnuid  as  indicated  on  Fiijure  21: 

1.  Obtain  the  elastic  solution  t.o  a  luniK’]  loaded  witii 
tile  O.B  kuar  (far-tH'ld  uniaxial  strain)  stres.ses. 


Figure  21.  Successive  Steps  in  Quasi-Static  Numerical  Simulation  Technique. 


2.  Using  tne  elastic  solution  as  initial  conditions, 
activate  the  elastic-plastic  material  model  and 
calculate  the  dynamics  until  a  steady  state 
solution  is  oiatained. 

3.  Using  this  elastic-plastic  solution  as  initial 
conditions,  a  calculation  is  performed  with  the 
slide  lines  introduced  to  model  the  joints  and 
permit  block  motion  and  deformation.  Also,  the 
material  may  fail  in  tension  if  the  tensile 
stresses  exceed  500  psi  in  any  principal 
direction.  The  dynamic  calculation  is  continued 
until  tensile  failure,  block  motion  and  tunnel 
deformation  have  stopped,  and  a  steady-state 
solution  is  thereby  obtained. 

4.3  JOINT  SLIDING 

Peak  slide  displacements  of  I'i  of  the  tunnel  diameter  are 
predicted  in  tlie  quasi-static  solution.  In  a  full  sized  tunnel, 
with  a  diameter  of  10  feet,  this  corresponds  to  1.2  inches  of 
sliding  . 

Figure  22  sliows  the  slide  displacement  along  the  joint  as  a 
function  of  distance  (scaled  by  tunnel  radius)  from  the  tunnel 
axis.  Tlie  slitiintj  is  a  maximum  at  about  one  radius  along  the 
joints  nearest  the  tunnc'i  crown  and  invert  (Joints  3  and  9  in 
Figure  22)  . 

4.4  TLNSILF,  CRACKINti 

The  tensile  cracking  is  closely  related  to  the  relative  slide 
d  i  sjj  lacements  across  the  joints.  Wlien  slipping  occurs  across  a 
joii\t,  ttie  plates  of  rock  separated  by  the  joints  act  as  beams. 
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The  presence  of  the  unlined  tunnel  removes  localized  support, 
and  beam  bending  leads  to  tensile  stresses  and  fracture. 

The  crack  pattern  obtained  in  the  numerical  solution  is 
shown  in  Figure  23.  The  "beams"  terminating  in  the  opening 
show  significant  cracking  of  the  unsupported  cantilevered  ends. 

4.5  TUNNEL  CLOSURE  AND  COMPARISON  WITH  EXPERIMENTAL  DATA 

The  tunnel  closure  predictions  are  shown  on  Figure  24.  The 
closure  along  the  vertical  axis  (i.e.,  crown-invert)  is  approxi¬ 
mately  0.5i.  The  maximum  closure  of  1.2%  is  located  at  0  =  30°. 

Figure  24  shows  a  comparison  of  the  actual  and  predicted 
tunnel  closure  (as  a  %  of  the  diameter)  versus  tunnel  angle.  The 
agreement  is  ciuite  good  except  near  the  crown-invert  (0  =  0°  and 
180°).  The  experimental  data  show  a  1.6%  closure  at  the  crown- 
invert  compared  to  a  0.5%  predicted  value.  This  difference  is 
probably  due  to  the  moderate  sliding  on  the  joints  and  associated 
tensile  failure  and  relaxation  which  occurs  in  the  calculations, 
but  not  in  the  experiment. 

In  the  ciuasi-stat ic  calculation,  tensile  cracking  occurs 
after  the  joints  arc  "unlocked"  and  joint  sliding  begins  (see 
Section  4.4  and  Figure  22).  Tensile  failure  is  not  evident  in 
the  post-shot  experimental  model  (I’igure  25).  If  the  metal  casing 
(wliich  is  not  modeled  in  the  calculation)  reduced  the  sliding  on 
the  joints,  or  if  the  actual  joints  arc  stronger  (higher  cohesion 
and/or  higher  coefficient  of  friction)  than  used  in  the  calcula¬ 
tion,  then  there  would  be  less  block  motion  in  the  experiment, 
and  therefore  there-  would  also  be  less  tensile  failure. 

Figuire  25  sh(jws  the  post-shot  experimental  model.  Note  tlie 
de  forrt’a  t  ion  of  th<’  tliin  (.003  inch)  metal  casing  (or  liner) 
surroundinu  thc'  frxim.  Tliis  casing  ii'ay  liave  inhibited  sliding 
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vl  ion  ;  t  h>-‘  joints.  .  ^  :  j'o  i  ■  K  ; 

■  •,  1  -  Tlius ,  "block  motion"  has 

occurri.'ii  as  f- rovlictovi .  A  quantitative  uescri[jtion  from  the 
snnii  i- sea  I  e  luo^iei  ol  tfie  amount  of  relative  dis[jlaccment  across 
tile  joints  is  niijlily  desirable  for  comparisons  v;ith  the  WAViJ-L 
code  predictions  of  I’lquro  2d. 
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APPENDIX  A 

JOINT  ACTIVATION  CONDITIONS  NEAR  A  CYLINDRICAL  TUNNEL 


The  relationships  between  the  principal  stress  components 
('■'if  Oj),  joint  and  tunnel  angles  (t,  i:i),  and  constitutive  acti¬ 
vation  properties  will  be  obtained  using  a  Mohr  diagram  (Ref.  8). 
The  joint  constitutive  model  is  specified  by 


max 


i  + 

o 


po 


n 


or 


where 


max 


+  ( tan  +■ )  c  ^ 


(A-1) 


=  the  maximum  allowable  shear  stress  on  the  joint 

=  the  normal  stress  on  the  joint 
n 

=  the  joint  cohesion 
o 

=  tan  ;  is  the  coefficient  of  friction  and  ;  is 
the  angle  of  friction. 


Figure  A-1  (page  57)  is  a  Mohr  diagram  which  shows  (1)  the 
shear  stress  and  normal  stress  ■'  for  the  joint  constitutive 
model  (Equation  A-1),  and  (2)  the  possible  shear  and  normal  stress 
states  corresponding  to  the  principal  stresses  i  and  .  .  The 
following  sketch  (page  54)  indicates  the  geometrical  interpretation 
of  the  Mohr  diagram  in  terms  of  the  stresses  on  a  joint  near  a 
tunnel  surface: 
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Doing  tho  algebra  and  solving  for  Op,  we  find, 


/ sin(2v'-;)  +  sin  ;a 

‘ \ sin ( 2 V- i )  -  sin  v 


o  / _  sin  1 _ 

li  lsin(2v-v)  -  sin 


Equation  A- 3  relates  the  two  principal  stresses  such  that 
an  oriented  joint  (with  angle  v  between  the  tunnel  normal  and 
the  ]oint)  will  lead  to  joint  activation.  A  more  physical  and 
useful  relationship  for  a  plane  stress  wave  engulfing  a  tunnel 
can  be  obtained  on  a  ,r  ,  rf  i--,- . 

On  an  tunnel  surface,  one  principal  stress  component 

vanishes;  i.e.,  ci  =0.  Thus,  the  tangential  stress,  ^2  r  on  the 
tunnel  surface  for  joint  activation  is 


\  2 1  cos  ; 

_ sin  : _ _  _ o _ I _ 

s  i  n  ( 2  ,  -  ;  )  -  sin  ~  sin(2v-;)  -  sin 


(A-4) 


When  a  plane  wave  engulfs  a  tunnel,  however,  the  tangential  sc 
varies  witn  tlie  angle  on  the  tunnel  surface  as  illustrated  in  the 
following  sketch: 


=  ^(■■)  ■ 


[ 1  +  Iros  (i  -  )  1 


Kf  1  at  i  on.sh  i  [)  l!et  ween  iaiifU'titial  Stress,  ■  , 

and  toadin;;  Stress,  a  on  luniu' 1  Surlaee 

I  oad 


In  the  sketch,  A(3)  is  the  amplification  factor  of  the 
Ofoad  '^3ve .  A(3)  depends  on  the  angle  3  on  the  tunnel  surface 

and  on  the  elastic  Poisson's  Ratio  (v)  of  the  media.  The  deri¬ 
vation  for  A(B)  can  be  obtained  by  superposition  of  elastic 
solutions  derived  in  Reference  8.  Figure  5  (in  body  of  report) 
shows  A(3)  for  v  =  .25  and  .4. 

Since  0  i  =0  and  Oo  =  A(3)  ,  Equation  A-4  shows  that 


load 


2t  cos  (1) 

_  0  7  _  _ o _ _ _ 

A(3)  [sin(2i);-^)  -  sin  4)]A(3) 


where 


A(3)  =  2 


)  1  +  2cos (23)1 


i-2v 


This  is  the  activation  stress  load  relationship  (Equation  1) 
used  in  the  report. 


(compress i ve ) 
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